Metal Surface Energy: Persistent Cancellation of Short-Range 
Correlation Effects beyond the Random-Phase Approximation 
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The role that non-local short-range correlation plays at metal surfaces is investigated by analyzing 
the correlation surface energy into contributions from dynamical density fluctuations of various two- 
dimensional wave vectors. Although short-range correlation is known to yield considerable correction 
to the ground-state energy of both uniform and non-uniform systems, short-range correlation effects 
on intermediate and short-wavelength contributions to the surface formation energy are found to 
compensate one another. As a result, our calculated surface energies, which are based on a non-local 
exchange-correlation kernel that provides accurate total energies of a uniform electron gas, are found 
to be very close to those obtained in the random-phase approximation and support the conclusion 
that the error introduced by the local-density approximation is small. 

PACS numbers: 71.15.Mb, 71.45.Gm 
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I. INTRODUCTION 

The widely-used Kohn-Sham formulation of density- 
functional theory (DFT)EI requires approximations to the 
exchange-correlation (xc) energy E xc [n(r)]. The simplest 
approximation to this functional is the so-called local- 
density approximation (LDA), where E xc [n(r)] is given 
at each point by the xc energy of a uniform electron gas 
at the local density. This approximation was found to 
be remarkably accurate in some rather inhomogeneous 
situations ,□ and its widespread use in condensed-matter 
physics led to the early success of DFT. 

Hence, it is important that the LDA be tested against 
benchmark systems, such as the jellium surface, and that 
new functional be developed. Nevertheless, more than 
30 years after Lang and Kohn reported the first sel: 
consistent LDA calculation of the jellium surface energy,l 
the question of the impact of non-local xc effects on the 
surface energy and their interplay with the strong chacge 
inhomogeneity at the surface has remained a puzzlefl'tJ 
The simple LDA and more advanced density functionals 
such as generalized gradient approximations (GGA's)3 
and meta-GGA'sQ all predict the same jellium xc surface 
energy within a few percent, but show no such agreement 
with the available wave-function based methods: Fermi 
hypernetted chain (FHNC)eI and diffusion Monte Carlo 
(DMC)Ja see Table I of Ref. [jlj. 

An alternative formally exact way to find the xc en- 
ergy of an arbitrary inhomogeneous system is provided 
by the adiabatic connection formula and the fluctuation- 
dissipation theorem.EU Within this approach, the ex- 
change energy is fully determined from the exact Kohn- 
Sham (KS) orbitals and the correlation energy is ob- 
tained in terms of the xc kernel J xc Ji3 In the random- 
phase approximation (RPA), f xc is taken to be zero. Full 
RPA or corrected-REA calculations are now feasible not 
only for bulk jelliumt 3 ! but also for jellium surfacesli-3 and 
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In this paper, we take a non-local xc kernel that pro- 
vides-accurate ground-state energies of a uniform electron 
gaSjE 3 ! and evaluate the jellium surface energy through 
the use of the adiabatic connection formula. We analyze 
the correlation surface energy into contributions from 
dynamic density fluctuations of various two-dimensional 
wave vectors, and find that short-range xc effects on 
intermediate and short-wavelength fluctuations nearly 
compensate. Hence, while RPA is known to be a poor 
approximation for the total correlation energy, our calcu- 
lations show that it is a surprisingly good approximation 
for those changes in the correlation energy that arise in 
surface formation. This js in contrast with FHNC and 
DMC slab calculations pu which predict surface energies 
that are significantly higher than those obtained either 
in the LDA0 or in a fully non-local RPAH 

The FHNC variational equations have been shown to 
provide, in the homogeneous limit, reasonable agreement 
with the known properties of a uniform electron gas,E3 
and DMC calculations are often regarded as essentially 
exact.113 Furthermore, one may expect that when applied 
to non-uniform systems these wave-function-based ap- 
proaches will lead to results whose accuracy is compa- 
rable to the high accuracy obtained for uniform systems. 
Nevertheless, we show that surface formation energies 
obtained from slab calculations either by a linear fit in 
the slab thicknessEa or as differences between slabuener- 
gies and an independently determined bulk energyO may 
result in substantial imprecision, and conclude that wave- 
function-based estimates need to be reconsidered. 



II. THEORETICAL FRAMEWORK 

We consider a jellium slab normal to the z axis, which 
is translationally invariant in the surface plane. Subtract- 
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ing from the slab energy the corresponding energy of a 
uniform electron gas and using the adiabatic connection 
formula, one obtains the xc surface energy 



d(q/k F )^ c 



(1) 



q being a wave vector parallel to the surface, kp is the 
Fermi momentum and 7^° represents the surface energy 
associated with the good quantum number q: 

lg C = / ' dz J dz' n(z) v q (\z - z'\) 

x / d\[nl%{ Z ,z')-nl%{\z-z'\)], (2) 
Jo 



v g (\z - z'\) = (2ne 2 /q) exp (-q\z 



being the Fourier 



transform of the bare Coulomb interaction. n xc x {z,z') 
and 



q,X\ 

represent Fourier components of the 
xc-holc density of a fictitious jcllium slab at coupling 
strength Ae 2 and the corresponding xc-hole density of a 
uniform electron gas of density n — kp/3ir 2 , respectively. 
In the LDA, the xc surface energy is obtained by simply 
replacing n xc x {z,z') by the xc-hole density of a uniform 
electron gas of density n(z). A parametrization of the 
uniform-gas-.xc-hole density has been given by Perdew 
and Wang,Ei which yields-the DMC ground-state energy 
of a uniform electron gasO 

According to the fluctuation-dissipation theorem, 

H f°° 

n q%( z > z ') = rr / dwx q ,\{z,z'\iu))- 8(z- z'), 

(3) 

where Xq,\( z , z ';u>) is the interacting density-response 
function. ' Time-dependent DFT (TDDFT) rfj^ws that 
this function obeys the Dyson-type equationE3t3 

Xq,\(z,z';ui) = x°(z,z';uj) + J dzi J dz 2 Xq(z,zv,u)) 
x [Xv q (\zi - z 2 \) + f^ x [n]( Zl ,z 2 ;oj)] Xq,\{z2,z';w), (4) 

where Xq( z , z'; u) is the density-response function of non- 
interacting KS electronscj and / ? ™[ri](z, z'; lo) involves 



the functional derivative of the KS xc potential at cou- 
pling constant A. 

Using the coordinate-scaling relation for the A- 
dependence of the xc kernel derived in Ref. O, we find 

f^ x [n(z)](z,z';w) = f*f x [A -3 n(z/ A)] (Xz, Xz'; lo/X 2 ), 

(5) 

where f xc [n](z, z'; lo) is the xc kernel at A = 1. In order 
to derive an approximation for this quantity, we assume 
that the density variation [n(z) — n(z')] is_small within 
the short range of f x c [n] (z, z'\ lo) and writell 3 ] 

/*>] (z, z'; w) = f* c ([n(z) + n{z')] /2; \z - z'\;uj), (6) 

where fq C (n; \z — z'\;u>) is the Fourier transform of the 
xc kernel f xc (n; k, lo) of a uniform electron gas of density 
n. Here k = (q, k z ) represents a three-dimensional wave 
vector. 

III. RESULTS AND DISCUSSION 



We have carried out simplified surface-energy calcu- 
lations with f xc (n; \z — z'\;co) replaced by f xc (n;k — 
q,Lo)S(z — z') [thus assuming that the dynamic den- 
sity fluctuation is slowly varying in the direction per- 
pendicular to the surface] and using _the parametriza- 
tion of Richardson and Ashcroft for f xc (n; k, lo)^3 and 
have found that neglect of the frequency dependence of 
the xc kernel does not introduce significant errors. We 
have also carried out adiabatic LDA (ALDA) surface- 
energy calculations with f xc (n;\z — z'\;lo) replaced by 
f xc (n;k = 0,lo = 0)5(z — z') [thus assuming that the 
dynamic density fluctuation is slowly varying in all di- 
rections] , and have found that the spacial range of the xc 
kernel cannot be neglected. .These conclusions also apply 
to the uniform electron gas.E3 

Hence, we neglect the frequency dependence of the 
xc kernel and exploit the accurate DMC calculations re- 
ported in Ref. [27| for the static xc kernel of a uniform 
electron gas. A parametrization of this data satisfying 
the known small- and lacae-wavelength asymptotic be- 
havior has been reported,E3 which allows us to write 
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where C, B, g, a, and (3 are dimensionless functions of 
the electron density (see Ref. ^8|) , n — kp/3n 2 , and z — 
z — z'. The finite q — > limit of Eq. (Q) will be dominated 
by the q — > divergence of v q (\z— z'\), making RPA exact 
in this limit. In the large-g limit, where short- wavelength 



excitations tend to be insensitive to the electron-density 
inhomogeneity, introduction of Eq. (Q) into Eq. (||) is 
expected to yield an xc kernel that is essentially exact. 

If the interacting density-response function 
X q ,\{z, z'\ lo) entering Eq. (||) is replaced by Xq( z , z']tv), 
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FIG. 1: Wave-vector analysis 7J of the correlation surface 
energy for a jellium slab of thickness a = 7.17r s and r s — 
2. Thin-solid, dashed, and dotted lines represent Perdew- 
Wang, RPA and Corradini-based LDA calculations, respec- 
tively. The thick-solid line and the open circles represent RPA 
and Corradini-based non-LDA calculations, respectively. Our 
'best' non-local calculation (open circles) provides the exact 
small-g and large-g limits. 

Eqs. (|]) and (||) yield the exact exchange surface energy, 
as obtained in Ref. [uj Here we focus our attention on 
the correlation surface energy, which for comparison we 
also calculate in LDA by replacing n c q A in Eq. (§) by the 
uniform-gas correlation-hole density at the local density 
n{z). 

Fig. 1 shows the wave-vector analysis 7^ of both LDA 
and non-LDA correlation surface enecgies of a jellium 
slab of thickness a = 7.17r s and r s = 2zB First of all, we 
focus on our LDA calculations, which have been carried 
out either by using the uniform-gas correlation-hole den- 
sity with f* c (n; \z - z'\) = (RPA-based LDA) or the 
xc kernel of Eq. (Q) (Corradini-based LDA), or by us- 
ing the Perdew-Wang (PW) parametrization of Ref. |2(]. 
We observe that in the long-wavelength limit (q — > 0) 
both RPA and Corradini calculations coincide with the 
PW parametrization. At shorter wavelengths, the Corra- 
dini scheme predicts a substantial correction to the RPA 
and accurately reproduces the PW wave-vector analysis 
of the correlation energy. Hence, armed with some con- 
fidence in the accuracy of our choice of the xc kernel, 
we apply it to the more realistic non-local scheme de- 
scribed above. Fig. 1 shows that our non-LDA beyond- 
RPA-(,'best' non-local) calculation coincides in the q — > 
limited with the non-LDA RPA, which is exact in this 
limit. In the large-g limit, local and non-local calcula- 
tions coincide, and our 'best' non-local calculation accu- 
rately reproduces the PW-based LDA (thin solid line), 
which is expected to be esentially exact in this limit. 

Consequently, our 'best' non-local calculation provides 
both the exact small-g limit, where LDA fails badly, and 
the exact large-g limit, where RPA is wrong. The LDA 
largely underestimates our non-local correlation surface 



TABLE I: Non-local xc (<j xc ) and total (a) surface energies 
and their local (LDA) counterparts, as obtained from Eqs. ([j]) 
and (^) with the non-local xc kernel of Eq. (^). Also shown 
are non-local RPA xc surface energies (orpa) and PW-LDA 
total surface energies (cpw— lda)- Tiny differences betwepi 
these RPA xc surface energies and those reported beforet-1 
are entirely due to differences in the parametrization of the 
xc potential. Units are erg/cm 2 . 
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energy, but Fig. 1 shows that the difference between RPA 
and beyond-RPA 7^ (the short-range part of the corre- 
lation surface energy) is fairly insensitive to whether the 
LDA is used or not. This supports the assumption made 
in Ref. [u] that the short-range part of the correlation en- 
ergy can be treated within LDA or GGA. However, short- 
range xc effects on intermediate and short-wavelength 
contributions to the surface energy tend to compensate, 
and this cancellation happens to be even more complete 
than expected from LDA or GGA. Thus, our non-local 
scheme yields surface energies which are still closer to 
RPA than is the RPA+ of Ref. |l} 

To extract the surface energy of a semi-infinite 
medium, we have considered three different values of the 
slab thickness: the threshold width at which the n = 5 
subband for the z motion is completely occupied and the 
two widths at which the n — 5 and n — 6 subbands are 
half occupied, and have followed the extrapolation proce- 
dure of Ref. In Table I we show our extrapolated local 
(LDA) and non-local surface energies, as obtained from 
Eqs. (0) and @ cither with |? c (n; \z - z'\) = (RPA) 
or with the xc kernel of Eq. (Q) . These calculations in- 
dicate that the introduction of a plausible non-local xc 
kernel yields short-range corrections to RPA surface en- 
ergies that are negligible. For comparison, also shown 
in Table I are PW-LDA surface energies, as obtained ei- 
ther from EqSj— ,(|l|) and (||) with the PW uniform-gas 
xc-hole densityE3 or from the PW parametrization of the 
uniform-gas xc energy.o Corradini and PW-based LDA 
surface energies (ctlda an d opw-lda) are found to be 
very close to each other, and we expect our Corradini- 
based non-local surface energies (a) to be close to the 
exact jellium surface energy, as well. 

We close this paper with an analysis of the avail- 
able wave-functioa-bascd surface-energy calculations. 
Krotscheck et al^B considered slabs of four different 
thickness a, and obtained both a bulk energy per par- 
ticle and a surface energy a from the FHNC slab 
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energy per particle e as a function of the particle number 
per unit area na, by a linear fit s(na) = + 2a /ha 
that becomes exact in the limit of infinite thickness. 
These authors showed that the extrapolated bulk en- 
ergies (Eoo) and those obtained from a separate FHNC 
bulk calculation (e) agree within about 1%, and claimed 
that this comparison lent credibility to their numerical 
treatment. However, these small differences in the bulk 
calculation yield an uncertainty in the surface energy 
Act = (e — £00)^0/2, which for r s = 2.07 and 4.96 can be 
as large as 280 and 11 erg/cm 2 , respectively. Moreover, 
since for the smallest /largest width under study and due 
to oscillatory quantum-size effects the quantity e(na) is 
larger/smaller than expected for a semi-infinite medium, 
the extrapolated bulk and surface energies are found to 
be too negative and too large, respectively. 

Li et alx3 calculated the fixed-node DMC surface en- 
ergy of a jellium slab with r s = 2.07 and found a = 
—465 erg/cm 2 , which is ~ 40 erg/cm 2 larger than the 
RPA value. They also performed LDA calculations with 
either the Wigner or the Ceperley-Alder form for the 
uniform-gas xc energy, and found LDA surface energies 
that are also about 40 erg/cm 2 larger than the corre- 
sponding LDA surface energies of a semi-infinite jellium, 
which suggests that finite-size corrections might bring 
the DMC surface energy into close agreement with RPA. 
These fixed-node DMC calculations were extended by 
Acioli and Ceperley to study jellium slabs at five dif- 
ferent densities,B but these authors extracted the sur- 
face, energy from release-node-,bulk energies. Both Li et 
alx3 and Acioli and Ceperleycl claimed that the release- 
node correction of the uniform electron gas at r s = 2.07 
is 0.0023 eV/electron, and argued that this correction 
would only yield a small error in the surface energy. How- 
ever, the unpublished uniform-gas fixed-node energy re- 
ported and used in Ref. [52] is actually 0.0123 eV/electron 
higher than its release-node counterpart; hence, by com- 
bining fixed-node slab and release-node bulk energies 
Acioli and Ceperley produced for r s = 2.07 a surface 
energy that is too large by 138 erg/cm 2 . Furthermore, 
had these authors used fixed-node bulk energies (see, e.g., 
Rcf . |33|) , they would have obtained surface energies that 
are close to RPA. 



IV. SUMMARY AND CONCLUSIONS 

We have investigated the role that short-range correla- 
tion plays at metal surfaces, on the basis of a wave- vector 
analysis of the correlation surface energy. Our non-local 
calculations, which are found to provide the exact small 
and large-? limits, indicate that a persistent cancellation 
of short-range correlation effects yields surface energies 
that are in excellent agreement with RPA, and support 
the conclusion that the error introduced by the LDA 
is small. Although this conclusion seems to be in con- 
trast with available wave-function based calculations, we 
have shown that a careful analysis of these data might 



bring them into close agreement with RPA. This is con- 
sistent with recent work, where jellium surface energies 
extracted from DMC calculations for jellium spheres are 
also found to be close to RPAo 

We have found that the RPA xc surface energy displays 
an error cancellation between short- and intermediate- 
range correlations. A different (and less complete) error 
cancellation between long- and intermediate-range xc ef- 
fects explainst2l why the LDA works for the surface en- 
ergy. The GCIA corrects only the intermediate-range 
contributions,!^ and so gives surface energies slightly 
lower and less accurate than those of LDA. Usually GGA 
works better than LDA, but not for the jellium surface 
energy where long-range effects are espj ecialhr .i irn portant. 
In the present work, as in four others ,E3cac3E3 we have 
found that the jellium xc surface energy is only a few 
percent higher than it is in LDA. These closely-agreeing 
methods include two different short-range corrections to 
RPA (present work using a non-local xc kernel and Ref. 
|3~i| using an additive GGA correction) , a long-range cor- 
rection to GGA (Ref. [Tcj), extraction of a surface energy 
from DMC energies for jellium spheres (Ref. [54]), and 
a meta-GGA density functional (Ref. [?]). The corre- 
sponding correction to. the LDA or GGA surface energy 
has been transferred^ successfully to the predictioa-of 
vacancy formation energiestH and works of adhesion £3 

We have almost reached a solution of the surface- 
energy puzale, but one piece still does not fit. Krotscheck 
and KohnEi examined a " collective RPA" which approx- 
imates our full RPA, and also used several xc kernels 
to correct for short-range effects. When they used an 
isotropic xc kernel derived from the uniform gas, in the 
spirit of our Eq. (6), they found surface energies, very 
close to RPA, as we do. When they used FHNCjS cor- 
responding to an anisotropic (but to the eye not very 
different) kernel constructed explicitly for the jellium sur- 
face, they found a large positive correction to the RPA 
surface energy, amounting at r s = 4 to as much as 35% 
of the RPA o xc or 60% of the RPA total a. The story 
of the jellium surface energy cannot reach an end until 
their work is reconciled with the other work on this sub- 
ject. The DMC energies of jellium slabs should also be 
re-considercd.c3 

Measured surface energies of real metals have been 
compared with calculations in Refs. f4l|,[42 43, However, 



experimental and calculational uncertainties and differ- 
ences between jellium and real-metal surfaces seem to 
preclude a solution to the surface-energy puzzle from 
these comparisons. 
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